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ABSTRACT 

We apply methods developed by Lovelace, Lipatov, and Kirschner to evalu- 
ate the leading Regge trajectories a{t) with the quantum numbers of nonexotic 
quark-antiquark mesons at A'^c = cxd in the limit t ^ —oo where renormalization 
group improved perturbation theory should be valid. We discuss the compatibil- 
ity of nonlinear trajectories with narrow resonance approximations. 



-k Work supported in part by the Department of Energy, contract DE-FG05-86ER-40272. 



It is unlikely that Quantum Chromodynamics, the consensus theory of strong 
interactions, can be exactly solved with realistic values for all parameters. How- 
ever, asymptotic freedom allows the application of weak coupling techniques such 
as perturbation theory to obtain the predictions of QCD for processes controlled 
by short distance dynamics. Besides high momentum transfer collision phenom- 
ena, one can hope to use such weak coupling techniques for computing the mass 
spectrum of hadrons containing only very heavy quarks. But for hadrons contain- 
ing light quarks and also for glueballs, strong coupling dynamics is unavoidable. 

It is reasonable to first confront these strong coupling issues making as many 
simplifying idealizations as possible. Thus we might idealize light quarks to mass- 
less quarks and delete the heavy quarks all together. Since m^, rud « ^qcd this 
idealized theory should give an excellent approximation to the dynamics of up 
and down quarks, which is to say all ordinary hadronic matter. These ideal- 
izations leave us with a theory with no further free parameters. Unfortunately, 
making the quarks massless does not simplify the dynamics enough for an ana- 
lytical treatment: the S'-matrix is non-trivial in all channels including those with 
particle production. Furthermore, the bound state spectrum includes all nuclei 
as well as the lowest mass hadron in each fiavor sector. 

That is why 't Hooft's idea of exploiting the AT^ — > oo limit '^'is so attrac- 
tive. In this limit the scattering amplitudes involving hadrons vanish and in 
lowest nonvanishing order are meromorphic in the channel invariants, just as the 
tree approximation to a quantum field theory. Nor do the nuclei bind in this 
limit. Thus an exact solution in this limit really would be significantly simpler 
than the exact solution at N^. = 3. Even if the definitive tests of QCD must 
come from large scale computation, a successful analytical understanding of all 
hadronic phenomena to 30% accuracy would be very desirable. Unfortunately, 
to date efforts to evaluate the large Nc limit have failed: at least with available 
methods the limiting theory seems almost as intractable as the finite Nc theory. 
Nonetheless, we think it is worthwhile to develop as much insight into the nature 
of the hoped for solution as possible. 
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String theory started as an effort to build exactly the sort of approximation 
to strong interaction dynamics that is provided by large Nc QCD. Since that 
approach led to the "wrong" answer, we should understand how the expected 
properties of QCD are different from those of string theory. In string theory the 
Regge trajectory functions astring{t) = Oi't + ckq where a' = 1/2ttTq with Tq the 
rest tension in the string, play a central role in the string scattering amplitudes: 
they appear directly as the arguments of the gamma functions which characterize 
the Veneziano four string function ^4(5, t) = ggfj.j^j^gr{—a{s))r{—a{t))/r{—a{s) — 
a{t)). The meromorphy of A4 in s and t follows directly from that of the gamma 
functions and the exact linearity of the Regge trajectories. Since the large 
hadron amplitudes are also meromorphic in s and t, the trajectory functions 
themselves should be good characteristics of the similarities and differences be- 
tween string theory and QCD. Also they might carry some hints about the 
solution of large Nc QCD. 

In this letter we study the Regge trajectories of large Nc QCD in the "me- 
son" channels {i.e. those interpolating the rotational states of quark-antiquark 

mesons), in the limit of large negative t where perturbative QCD should be ap- 
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plicable. We follow ideas and methods developed by Lovelace, Lipatov, and 

[4] 

Kirschner and Lipatov. These methods are essentially renormalization group 
improved calculations based on summing leading logarithmic contributions of 
Feynman graphs. Such methods can only give the trajectory functions in the 
weak coupling approximation. Since the coupling "runs" with the scale A(— t) = 
TVc^ff (— t)/47r^ ~ 12/11 \'n.{—t/ Aq(jjj) this means that we can obtain only the large 
negative t behavior of the Regge trajectories using these methods.* 

Our first task is to identify the leading logarithmic contributions to a scat- 
tering process involving the exchange of a qq pair. In any gauge theory the 
ladder diagrams, which iterate gauge boson exchange between two fermion lines, 

* Note that the large positive t behavior of the trajectories is characterized by the confining 
force and should be asymptotically linear aQcoit) ~ t/2TTk where kR is the confining term 
in the qq interaction energy. 
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contribute two powers of ln(s//i^) for each additional rung. Thus the leading 
logarithms are actually double logarithmic and dominate the single logarithms 
of renormalization: the leading log sums will therefore not include running cou- 
pling effects. Thus we proceed in two steps. First we evaluate the amplitudes to 
double logarithmic accuracy and then incorporate renormalization effects which 
make the coupling run in the second step. As shown in 1967 for QED, the first 
step typically leads to a fixed square root branch point in the angular momentum 
plane. For QED processes involving the exchange of total zero charge, the lead- 
ing double logs come only from the ladder sum, which produces a branch point 
located at J = ^Jlaji^ where ol ~ 1/137 is the fine structure constant. When 
non zero charge is exchanged, there are additional double logarithmic contributi- 
ions coming from soft "bremstrahlung" photons which either form crossed rungs 
in the ladder or Sudakov vertex corrections. Similarly, in QCD^ as discussed by 
Kirschner and Lipatov, the double logarithmic contributions of the basic ladder 
diagrams of Fig. 1 are supplemented by soft gluon bremstrahlung graphs. For- 
tunately for us, these additional diagrams are nonleading in the 1/Nc expansion 
when the ladder structure is "hooked on" to color singlet hadron vertices^ and 
so do not enter into our calculation of the Regge trajectories of large Nc QCD. 
Thus the double logarithmic sum is identical to the QED zero charge exchange 
case with A = Ncgl/^'n^ substituted for 2a /tt. 

Since renormalization effects are neglected in the leading double logarithmic 
approximation, the results depend on a fixed coupling constant A. Asymptotic 
freedom must at least make the location of the singularities in the angular momen- 
tum plane vary with t according to the replacement A A(— t), but actually the 
cut is changed into a distribution of Regge poles accumulating at as t — * — oo. 
This phenomenon was first uncovered by Lovelace for the case of (f)^ theory in 6 
space-time dimensions where the accumulation point is at J = —1. He analysed 
the Bethe-Salpeter (B-S) equation with a kernel improved to include the effects 

t Note that these hadron "form factors," represented by the left and right hand parts of Fig. 
1, involve on-shell mesons and cannot be calculated perturbatively. 



4 



of asymptotic freedom. This equation produced partial wave amplitudes with 
only pole singularities in the angular momentum plane. Since he only consid- 
ered the case t = 0, his results for the pole locations (Regge intercepts) were 
untrustworthy. (The low momentum theory is a strong coupling one and the B-S 
equation is not valid there.) This shortcoming was removed by Kirschner and 
Lipatovf' who incorporated t dependence in leading order and obtained a{t) for 
large negative t instead of a{0). With large enough t, the effective coupling is 

[3] 

weak, justifying the B-S equation. Earlier, Lipatov had obtained similar results 
for the asymptotic behavior of the pomeron (glueball) trajectory in QCD. In 
this note we find the corresponding asymptotic behavior of the qq trajectories in 
large QCD. 

In order to incorporate the running coupling, we consider the B-S equation 
which sums the ladder subgraphs in Fig. 1. We represent the Green's function for 
the ladder subgraphs as a matrix "ifab in the Dirac indices of the qq lines coming 
in at the left. We apply the Dirac operators to these two lines in the coordinate 
representation to obtain (for simplicity we take all niq = 0) 

7 • di'^ixi, X2; yi, 2/2)7 • ^2 = -S{xi - yi)S{x2 - ^2) + K^T2)Y'^YdiJ.Axi2) 

where we have followed Lovelace's treatment of (pQ, replacing the coupling con- 
stant by the running coupling X{x~'^) ~ —12/11 Ih^x^Aq^j;)). In a general covari- 
ant gauge we define the coordinate space propagator by d^i,{x) / Air"^ with 



The B-S equation is not gauge invariant, but violations of gauge invariance will 
be small for weak coupling. Thus if we only keep leading order answers, our 
results should be gauge invariant. We keep C, arbitrary so we can confirm this. 




+ 



(1 - Ox^,x„ 
x^ 



We expect this equation to be accurate when A << 1 i.e. for x\2A%cd 



« 1. The 
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singularities in the t channel are controlled by the solutions of the homogeneous 
equation 

It is convenient to work with CM and relative coordinates r = {xi + X2)/2 and 
p = Xi — X2 and to Fourier transform with respect to r, whose conjugate variable 
is q so that = —t. Then the homogeneous equation reads 

We can only make use of this equation for small p and large q when the effective 
coupling associated with both scales in the problem is small. 

For qp << 1, but p~^ » ^qcd, the B-S equation becomes quite manage- 
able. This limit reduces it to the g = case, with its 0(4) symmetry. As in 
Ref. 4 one can conveniently consider the partial waves in this limit by making 
the ansatz 



(p-0 



£-1 



P' 



/(|p|)+e-7^(|p| 



(1) 



ipi 

with a fixed light-like four- vector, so that one is forming traceless symmetric 
tensors of rank I. Plugging this ansatz into the B-S equation then yields the pair 
of equations 

[{V - £)2 _ 4]/ + 2[D -i\{V-i- 2\g =X{p-')[Cf - (1 - C)^] 

2 

2£[D + l]f-[V-£]'g ^X{p~')g. 

Here V = \p\d/d\p\ — d/dlii{\p\A) is simply the scaling derivative in the magni- 
tude of p. 

Eq.( 2) is a spinor version of the small p B-S equation for the 0g theory 
analyzed in Ref. 2,4. With A(p~^) approximated by — 6/11 ln(|p|A), the Laplace 
transforms of A/, A^r satisfy a pair of first order differential equations in D = — 2zi>. 
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Unlike the single equation in the cj)\ case which can immediately be integrated, 
this pair of equations is equivalent to a single component second order equation 
which cannot be so readily solved. However, the leading angular momentum sin- 
gularity (with largest Re t) is controlled (for small A) by I and v small compared 
to unity. This is made clear by considering the determinant of the coefficient 
matrix on the l.h.s. of these equations: 

In the small p, v limit, the second of ( 2) shows that 2lf ^ (A(p^^) + ["D — P\^)g- 
Inserting this approximate form for / into the first equation and making the same 
approximations there gives {—T>'^-\-£'^)g fa X{p~'^)g the Laplace transform of which 
gives a first order equation in u. Thus we see that this leading Regge singularity 
is controlled by equations independent of the gauge (, as we anticipated for weak 
coupling. Looking back to the full gauge dependent small p equations, we notice 
that in Landau gauge {( = 0) one can eliminate / in favor of ^ in a p independent 
way. Thus in this gauge the Laplace transformed equations are first order in v 
and can be directly integrated. 

We exploit this simplification by setting ^ = in the following. Solving the 
second equation for / and substituting in the first we obtain 

[f _ p2| ^ - V^]g =[4:-V^-f- 2e]X{p-^)g. (3) 

The equation for g is now quite similar to the (/>g case and we can repeat the 
steps in Ref. 4 to derive the asymptotic behavior of the Regge trajectories. The 
equation for g is first solved by writing it in the form 

where we have put V = —2iv. Since [i', — \n.{p^Pi?)] = —i we can interpret 
R = — ln(p2A2) — > idjdv^ integrate the equation and then transform back to the 
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coordinate representation to obtain 



dv exp < ivR — 



m 



4 + 4z/' - ^2 _ 2£ 

[£2 + 4zy'2][(£ + 2)2 + 4z/'' 



(4) 



This solution of the small p B-S equation is, in fact, a solution of the full B-S 
equation for q = 0. But of course the B-S equation is only a good approximation 
for large q and small p. In Ref. 4 the analogous solution for the 0g theory is used 
to gain information about the large q behavior of the trajectories by noting that 
for pq << 1 but p not too small one can have A(p~^) ~ for a large range 

of p (essentially because the scale dependence of A is only logarithmic). Thus 
instead of requiring regularity of the solution at i? = as in Ref. 2, the solution 
is matched to that of the B-S equation with a p independent coupling taken to be 
X{q^) « 1. For consistency of the weak coupling approximation this matching 
must be imposed at large R (small p). 

The large R behavior of ( 4) is exponentially damped for extremely large R 
but there is oscillatory behavior for R not too large. This can be extracted by 
finding the saddle points 



22/)' 




1 + e 



22R 



4 



+ 



3(4-£2_2£) e^{e + 2y 



UR 



16 

(5) 

which is 

consistent with large R provided \£\ « 1. In this regime, is well approximated 
by the saddle point evaluation 



We see that there are two real values of provided i? < 



£2(^+2)2 




cos 



- + uoR- — ^(z^o) 



where is the positive square root of ( 5) taken with the plus sign, and 



(6) 



$(z/o) 



2-£2-£ 1 
-tin ■ 



2iuo/£ 



+ 



£(£ + 3) 



-iln 



l-2iiyo/{£ + 2) 



8£{£ + l) l + 2iuo/£ 8{£ +!){£ + 2) l + 2iuo/{£ + 2) 
For sufficiently large R, the R dependence of vq can be neglected and ( 6) can 
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be matched to the small p solution of the B-S equation with constant coupling 
A(g^). Scale invariance of the finite q B-S equation implies that the solution is a 
function of p\q\. In the small p limit it therefore has the behavior 



2|„|2 



+ i^iln|p| |g 



(7) 



where vi is with 12/1 li? replaced by X{q^)- 

Replacing uq by ui in ( 6) and comparing to ( 7) gives the matching condition 

1 



\{q')<^{u^)-v^ 



5 TT 



(8) 



where r = integer. We shall find that the consistent small coupling solution of 
these equations gives i = 0{X^^'^) and Vi = 0(A^/^). Thus I'l/P = 0{X^^^^) and 
i/^/i'^ — 0(1). Neglecting all terms that vanish at zero A in the matching 

condition simplifies to 



S TT 



With these approximations A fa ^^(1 + Auf/l"^) so the leading term on the l.h.s. 
cancels and we are left with 



3P ~ ln(QVA2) V ^ 2 ^ 4 



Replacing P by X{q^), we thus find the Regge trajectory asymptotics 



aAR) 1-2 



11 
32 



-1 2/3 



(9) 



Notice that we have an infinite number of trajectories accumulating at in the 
limit t — > — oo, with almost identical behavior to those of the 0g theory obtained 
in Ref. 4. 
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The phase 6 is not determined from the smaU p dynamics considered so 
far. It must be determined by the dynamics at pq = 0{1). However, except in 
exceptional cases, 5 = tt in the hmit we are considering. This is because this 
limit involves fa in ( 7) . For u = the two behaviors p^^*'' are replaced by 1 
and In p. Generically, both behaviors will be present, and unless the coeflBcient 
of In p exactly vanishes, the behavior for u slightly different from zero must be 

N {^[{qp)-'''' - {qpfn + C2{qp)-'"' + C3(gp)+^' j 

with Ca finite at = 0. If the Inp term is present at = 0, then ci 7^ 
there and 5 = n at u = 0. This is analogous to the generic vanishing of phase 
shifts at zero energy. In that analogy the case ci = corresponds to a "zero 

energy resonance." For the 0^ case, 5 = tt was shown by explicit solution of the 

[4] 

constant coupling B-S equation using conformal invariance. In gauge theories, a 
conformal transformation changes the gauge condition, so the B-S equation, being 
gauge noninvariant, is scale invariant but not conformally covariant. Lacking an 
explicit solution, we can only state that it is likely, but not proven, that S = n 
for large A^^ QCD. 

We close with some comments about the significance of nonlinear Regge tra- 
jectories for large Nc QCD. There is a common belief'^ that narrow resonance 
approximations require exactly linear (or at worst polynomial) Regge trajecto- 
ries. However this conclusion depends on a maximal analyticity assumption that 
the trajectory functions are free of singularities in the t plane cut on the right 
at threshold branch points. Since Im the discontinuity across the thresh- 
old cut, is proportional to the resonance widths, the trajectories would then be 
entire functions in the limit of zero width resonances. We have seen that the 
qq trajectories of large Nc QCD approach constants as t — > —00, and confine- 
ment together with infinite Nc implies linear behavior as t — > -|-oo as well as no 
threshold branch points. The inescapable conclusion is that the maximal ana- 
lyticity assumption fails for the Regge trajectories of infinite Nc QCD and there 
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are additional singularities in the t plane. This is probably also true at A'^c = 3 
since there is no good physical basis for the absence of additional singularities. In 
Ref. 9 (for earlier models see also Ref. 10,11) one of us discussed some examples 
of narrow resonance models with nonlinear trajectories with algebraic branch 
points in the complex t plane. Some of these models have trajectories which are 
asymptotically linear for large positive t and approach constants at large negative 
t. Unfortunately, they approach these constants as an inverse power of t rather 
than an inverse power of ln(— t), so they are not candidates for large Nc QCD. 
Nonetheless, they do show that nonlinear trajectories are compatible with narrow 
resonances, and indicate a direction toward solving large Nc QCD. 
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Figure Caption A typical large Nc diagram contributing to meson scattering 
with the exchange of a qq ladder structure. The leading log approximation as 
s — > oo is the sum of graphs with an arbitrary number of gluon rungs represented 
by the vertical double lines. 



12 



